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Abstract 

We consider the performance of non-optimal hedging strategies in exponential 
Levy models. Given that both the payoff of the contingent claim and the hedging 
strategy admit suitable integral representations, we use the Laplace transform approach 
of Hubalek et al. iTTTi to derive semi-explicit formulas for the resulting mean squared 
hedging error in terms of the cumulant generating function of the underlying Levy pro- 
cess. In two numerical examples, we apply these results to compare the efficiency of 
the Black-Scholes hedge and the model delta to the mean- variance optimal hedge in a 
normal inverse Gaussian and a diffusion-extended CGMY Levy model. 

Keywords: Laplace transform approach, mean-variance hedging, delta hedging, 
Levy processes, model misspecification 

1 Introduction 

A basic problem in Mathematical Finance is how the issuer of an option can hedge the 
resulting exposure by trading in the underlying. In complete markets, the risk can be off- 
set completely by purchasing the replicating portfolio. In incomplete markets, however, 
additional criteria are necessary to determine reasonable hedging strategies. 
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A popular approach studied intensively in the literature over the last two decades is 
mean-variance hedging. Expressed in discounted terms, the idea is to minimize the mean 
squared hedging error 

( 7 r T \ 2 \ 

(1.1) 

over all in some sense admissible trading strategies Here, the random variable H is the 
payoff of the option, c is the initial endowment of the investor^ and S is the price process 
of the underlying. Since the stochastic integral ' $tdS t represents the cumulated gains 
from trading ( |1.1[ ) amounts to comparing the terminal value c + J Q T $ t dS t of the hedging 
portfolio and the option's payoff in a mean-square sense. Comprehensive overviews on the 
topic can be found in [fT9l l23l . For more recent publications, the reader is referred to Q 
and the references therein. In particular, semi-explicit representations of the minimal mean 
squared hedging error have been obtained in J6l [HI for exponential Levy models and in 
[TT4l [T5ll for affine stochastic volatility models by making use of an integral representation 
of the option under consideration. 

However, in practice delta hedging is still prevalent, where the hedge ratio is given by 
the derivative of the option price with respect to the underlying. Therefore, it seems desir- 
able to compute also the mean-squared hedging error of such a non-optimal strategy. For 
exponential Levy models, this has been done in the unpublished thesis IfTOll for continuous 
time under the restriction that the asset price process is a martingale and, more recently, by 
O in a discrete time setup, both using the approach of [fTTI . The contribution of the present 
study is the extension of the results of ifTOl to general exponential Levy processes and a 
larger class of hedging strategies. The resulting semi-explicit formulas for the correspond- 



ing hedging errors can be found in our main result, Theorem 4.2 With these results at hand, 
we also compare the performance of the Black-Scholes hedge and the model delta to the 
mean- variance optimal hedge in several illustrative examples. 

This article is organized as follows. Subsequently, we describe the setup for the price 
process of the underlying and for the payoff function. In Section [3} we introduce the class 



of hedging strategies to which our approach applies. These ^-strategies (cf. Definition 3.1 ) 
include in particular the Black-Scholes hedge and, more generally, delta hedges in expo- 
nential Levy models. In Section |4} we state the main theorem on the hedging error of 
A-strategies and provide a sketch of its proof on an intuitive basis. Finally, we illustrate this 
result by two numerical examples in Section [5} In the first, we compare the performance of 
the Black-Scholes hedge and the variance-optimal hedge in a normal inverse Gaussian Levy 
model whose parameters are inferred from a historical time series. In the second example, 
we study the hedging errors of the Black-Scholes strategy, the model delta hedge, and the 
variance-optimal strategy in a diffusion-extended CGMY model for parameters obtained 
by calibration to option prices. The technical proof of the main theorem is delegated to 



'Note that we suppose here that this quantity is given exogenously. Typically, it will be chosen to match 
some arbitrage-free price the investor receives for selling the option. However, it can also include other initial 
holdings or liabilities. 
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Appendix [A} 

For stochastic background and terminology, we refer to the monograph of Jacod and 
Shiryaev [fT3l . For a semimartingale X, we denote by L(X) the set of X-integrable pre- 
dictable processes and write ip • X for the stochastic integral of a process ip £ L(X) with 
respect to X. By (X, Y), we denote the predictable compensator of the quadratic covari- 
ation process [X, Y] of two semimartingales X and Y, provided that [X, Y] is a special 
semimartingale (cf. [12, comment after Theoreme 2.30]). 



2 Model and preliminaries 

In this section, we state our assumptions on the asset price process and the payoff. Note that 
we use the same setup as [fTTI for mean- variance hedging. 



2.1 Asset price process 

Let T > be a fixed time horizon, and denote by (fi, JP, (=^t) te r T ] , P) a filtered probabil- 
ity space. The discounted price process S of a non-dividend paying stock is assumed to be 
of the form 

S t = S e Xt , S eR +1 te[0,T], 
for a Levy process X. We demand that the stock price is square-integable, 

E(Sl) < oo, (2.2) 

which is a natural requirement when using the second moment of the hedging error as a risk 
criterion. The entire distribution of the Levy process X is already determined by the law 
of X\, which can be characterized in terms of the cumulant generating function k : D — )■ C, 
i.e., the unique continuous function satisfying 



E (e zXt ) 



e 



tK,(z) 



for z £ D := {z £ C : E(e Re( - z)Xl ) < oo} and t £ R+. Note that Condition (2.2) implies 

{z £ C : < Re (z) < 2} C D. 

Moreover, we exclude the degenerate case of deterministic S by demanding that X\ has 
non-zero variance, i.e., 

k(2) - 2«(1) ^ 0. 



2.2 Laplace transform approach 

In order to derive semi-explicit formulas in concrete models, we employ the Laplace trans- 
form approach, which is widely used in option pricing (cf., e.g., [|5l|2T]|) and applied by ifTTTl 
in the context of mean- variance hedging. The key assumption is the existence of an integral 
representation of the payoff function in the following sense. 
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Assumption 2.1. Let the payoff H of a contingent claim be of the form H = f (St) for 
some measurable function / : (0, oo) — > K, which admits the representation 

/■R+ioo 

/(s)= / s z p(z)^ (2.3) 

for p : C — > C and R £ R such that x >->■ + ix) is integrable and 

£ (e 2RXl ) < oo. (2.4) 



Note that Condition (2.4) implies H £ L (P), which is again a natural assumption in 



view of the problem at hand. 

Example 2.2. Most European options allow for an integral representation as in Assump- 



tion 



2.1 For example, the payoff function f(s) = (s — K) + of a European call with strike 



K > can be written as 

for arbitrary R > 1, cf. ifTTl Lemma 4.1]. In general, representations of this kind can be 
obtained by inverting the bilateral Laplace transform of the mapping x i-)- f(e x ) via the 
Bromwich inversion formula. More details and examples can be found in Q~H Section 4] 
and ED Chapter 3]. 

Henceforth, we consider a fixed contingent claim H = /(St) satisfying Assump- 
tion |2U 

In the following, we will represent several objects as integrals with respect to the weight 



function p from Assumption 2. 1 The following terminology allows to conveniently express 



that such integrals are well-defined. 

Notation 2.3. A measurable function h : (R + iR) — > C is called p-integrable if 



/oo 
\h(R + ix) | \p(R + ix) | dx < oo. 
■oo 



Analogously, a measurable function h : (R+iM.) x (i?+iR) — > C is called twice p-integrable 
if 



OO POO 



\h(R + ix,R + iy) \ \p(R + ix) \ \p(R + iy) \ dxdy < oo. 



— oo J — oo 



Moreover, for a more convenient notation we will always write Jn-i^ \h( z ) \ \p( z ) \ \dz\ for 
the integral \ h(R + ix) \ \p(R + ix)\ dx, and analogously for double integrals. 

Note that in particular, a bounded function is p-integrable since z ^ 1 is p-integrable 



by Assumption 2.1 
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3 A-strategies 



We now introduce the class of strategies for which we will compute the mean squared 
hedging error in Section |4} Moreover, we discuss the most prominent examples. 

As in , we focus on hedging strategies which allow for a similar integral representa- 



tion as the payoff function (cf. Section 2.2 ). 



Definition 3.1. A real-valued process ip is called A-strategy if it is of the form 

/■R+ioo 

(ft= ip(z) t p(z) dz 

J R—ioo 

with 

<p(z) t = S*: 1 g(z,t) 
for a measurable function g : (R + I'M) x [0, T] — > C such that 

1. 1 1 — y g(z, t) is continuous for all z E R + iR, 

f T 2 

2. z i — y I \g(z, s)\ 2 ds is p-integrable. 

Jo 

Remark 3.2. Condition [T] is required to ensure the existence of a unique solution of an 
ordinary differential equation (ODE) with right-hand side g(z, •). Condition [2j ascertains 
that ip is integrable with respect to S and that the cumulative gains J Q <ptdS t of the strategy 



possess a second moment (cf. Lemma A. 2). As a side remark, it also implies that any 



A-strategy is admissible in the sense of [|TTL Section 3]. 

To motivate this definition, we now show that these A-strategies generalize the so-called 
delta hedges, which are obtained by differentiating the option price with respect to the un- 
derlying in an exponential Levy model. We also recall that the optimizer of the mean- 
variance hedging problem ( |1.1[ ) is another special case if it is computed under a martingale 
measure for S. The most important concrete example for both is the Black-Scholes strategy, 
i.e., the hedge obtained by differentiating the Black-Scholes price. This strategy also allows 
to achieve perfect replication in the Black-Scholes model and hence minimizes < | 1 . 1 1 ) . How- 
ever, it is important to keep in mind that it leads to a non-trivial hedging error if applied in 
a different Levy model with jumps. 

Example 3.3. (Delta hedge) Computing the derivative of a price process with respect to 
the underlying in an exponential Levy model leads to a A-strategy. To see why this holds, 
denote by S t = S e Xt an exponential Levy process with driver X and associated cumulant 
generating function k under some martingale measure Q. Note that due to, e.g., model 
misspecification, the hedge may be derived in a model differing from the one where it is 
eventually applied, which is why we distinguish between the processes S and S. Using 
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Fubini's Theorem and the independence of the increments of X with respect to Q, the price 
process of the contingent claim with payoff function / in the model S is given by 

/ ~ \ / fR+ioo ^ \ r R+ioo 

E Q (f(S T ) & t ) =E Q [ S z T p(z) dz &t)= E Q (s z t JP t ) p(z) dz 

V J \ J R-ioo J J R-ioo V 7 



R+ioo 



S z e K(z)(T-t) p ( z j ^ 



R—ioo 



The delta hedge of the contingent claim in this pricing model is then given by ip A (S t ^,t) 
for 



a pR+ioo pR+ioo 

<p±(s, *):=!-/ s z e^ T -^ p(z) dz = / s^ze^^ p(z) dz, 

& S J R-ioo J R-ioo 

provided that the derivative exists and that integration and differentiation can be inter- 
changed. This is the case if z i-> ze K ^ T ~ t ^ is p-integrable for all t E [0, T), which is 
satisfied if either the the distribution of driver X is sufficiently regular or the payoff func- 
tion is smooth enough, cf., e.g., [8|. Using the resulting hedge in the model S yields the 
strategy 

/•R+ioo 

if A (S t .,t)= / St 1 ze^ T - t ^ V {z)dz, (3.5) 

J R—ioo 

which is a A-strategy if Condition [2] of Definition 3.1 is satisfied. 



Example 3.4. (Mean-variance optimal hedge in the martingale case) By the integral 
representation provided in [fTTT Theorem 3.1], the mean- variance optimal hedging strategy 
in an exponential Levy model is a A-strategy, provided that the corresponding asset price is 
a martingale. In this case, one can therefore use the results of the present paper to quantify 
the effect of model mis specification arising from using a mean-variance optimal hedge in 
another model. If the asset price process fails to be a martingale, the corresponding hedge 
contains a feedback term and therefore is not a A-strategy. Nevertheless, the results from 
the martingale case should typically serve as a good proxy, because numerical experiments 
using [11 , Theorems 3.1 and 3.2] supply compelling evidence that the effect of a moderate 
drift rate is rather small for mean-variance hedging. 

As stated above, the delta hedge and the mean-variance optimal strategy coincide in the 
Black-Scholes model. Moreover, the regularity conditions in Definition |3.1| are automati- 
cally satisfied in this case. 

Lemma 3.5. Let S be a geometric Brownian motion without drift, i.e., 

S t = S Q e-^ 2t +° w \ S eR+, te[0,T], 

for a constant a G (0, oo) and a standard Brownian motion W. Then the delta hedge and 
the mean-variance optimal hedge in the model S coincide and are given by the A-strategy 

/•R+ioo 

V f s = / ^-i^M-D(T-t) p{z) dz _ 

J R—ioo 



6 



Proof. We follow the lines of Example |3.3|and show that the necessary regularity condi- 



tions hold in the Black-Scholes case. First, the cumulant generating function of the driver 
of S is given by 

k(z) = -a 2 z(z — 1) 



for z G C. Using the arguments of Example |3.3[ we obtain the price process of H in the 
model S as 



E 



(f(S T ) ft) = / S z t e^ T ^p{z)dz. (3.6) 

Note that the conditional expectation exists and Fubini's Theorem can be applied because 
all exponential moments of the normal distribution are finite. For the further considerations, 
note that 

Re (k( z )) = l -a 2 (Re {zf - Re (z) - Im (zf) , 

and define 

g(z,t) := ^ (2)(T -* } 

for z G R + iR and t G [0, T]. Since the mapping x \-t xe~ ax ' 2 ', a > 0, is bounded on R + , 

z^\g(z,t)\< (e^'^Vl) (l^l + IIm^le-^ 21 " 1 ^ 2 ^) 

is bounded on R+iM. for fixed t G [0, T) and hence p-integrable. Consequently, the integral 
Jr^ 1 ™ s z ~ 1 g(z,t) p(z) dz is well-defined for s G (0, oo) and t G [0, T), and dominated 
convergence yields 



d_ 
ds 



R+ioo fR+ioo 

s z e n{z){T - t) p(z)dz= I s z ~ 1 g(z,t)p(z)dz. 

R—ioo J R—ioo 



Therefore the delta hedge of the price process in (3.6) is given by 



r> R+ioo 

tf:= I S z I 1 g(z,t)p(z)dz = rf s . 

J R—ioo 

Note that whereas this integral is not well-defined for t = T, the cumulated gains process 
y? A • St does not depend on the value of ^ . Since delta hedging leads to perfect replication 
in the Black-Scholes model, y? A is clearly mean-variance optimal. Let us now verify that 
(p A is indeed a A-strategy. Obviously, g(z, ■) is continuous for fixed z G R + iR. Moreover, 
we have 

£\g(z,a)\ 2 da< (e^"^ V l) (tR 2 + £ Im {zf e - 2 M*) 2 (T- s ) ds 

Elementary integration yields that the right-hand side is uniformly bounded for z G R + iR, 
which implies that the left-hand side is p-integrable. Thus Conditions [T] and [2] of Defini- 
tion |3J] are satisfied and we are done. □ 
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Remark 3.6. Note that ip BS is also a replicating strategy and in particular mean-variance 
optimal for geometric Brownian motion with drift, i.e, 



S t = S e^-^ 2 y +aWt for fi G (0, oo). 



Remark 3.7. It can be shown by basically the same arguments as in the proof of Lemma 3.5 



that the delta hedge in the sense of Example 3.3 exists, admits a representation as in (3.5 ) 



and is a A-strategy if the risk-neutral driver X has a Brownian component. 



4 Performance of a A-strategy 

We measure the performance of a strategy in terms of the resulting mean squared hedging 
error, i.e., the objective function used in mean-variance hedging. 

Definition 4.1. For a given initial endowment c G 1 and a strategy cp G L(S) such that 
(f • St G L 2 (P), we define the mean squared hedging error of the endowment/strategy pair 

(c, tp) as 

E((H-c-(f S T ) 2 ) ■ 

4.1 Main result 

In this section, we state the main result of this paper. For better readability, the proof is 
deferred to Appendix [A) 

Theorem 4.2. Consider an initial endowment c G R and a A-strategy <p of the form 

/■R+ioo 

(ft= ip(z) t p(z) dz 

J R—ioo 

with 

^z) t = St l g{z,t). 

Then the mean squared hedging error of the endowment/strategy pair (c, tp) is given by 

/•R+ioo pR+ioo 

E((H-c-tp- S T ) 2 ) = (w-c) 2 + / J(y, z) p(y)p(z) dydz, (4.7) 

J R—ioo J R—ioo 



where 



a(z, t) := ( 1 - k(1) e K ^ s - T) g(z, s) ds^j e^^, (4.8) 

"R+ioo 

w:= I S^a(z,0)p(z)dz, (4.9) 

R—ioo 
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and 



J(y, z) := («(y + z) - k{ V ) - k{z)) / S y +Z e^ +z > a(y , s)a(z, s) ds 

Jo 

- («(y + 1) - «(y) - «(1)) / S^V<>**a(y, s)^, a) 

</o 

- (k(z + 1) - «(z) - «(1)) / S y +z e K{lJ+z)s a(z,s)g(y,s)ds 

Jo 

+ (k(2) - 2«(1)) / S2f + V(» + *^(y, S )^,s)ds. (4.10) 
Jo 

Remark 4.3. The cumulant generating function k is often known explicitly, e.g., for nor- 
mal inverse Gaussian 0, variance gamma ifTTl and CGMY Hi processes or for the models 
introduced by Merton [|T8ll and Kou |[T6ll . Moreover, the time integrals in (|4.8|) and in (|4.10|) 



can typically be calculated in closed form, which means that the evaluation of the hedg- 
ing error ( |4.7[ ) usually amounts to numerical integration of a double integral with known 



integrand as in IfTTl Theorem 3.2] for the mean- variance optimal hedge. 



4.2 Sketch of the proof 



Besides giving a rigorous proof of Theorem 4.2 in AppendixfA] we present our approach on 



an intuitive level. The notation used in Appendix |A]is anticipated here, but it will be defined 
precisely in the corresponding places. 
To calculate 

E({H-c-ip'S T f), (4.11) 



we look for a martingale L with Lt = H — c — ip m St- Then we can rewrite (4. 1 1 ) as 



E((H-c-<p- S T ) 2 ) = E{Ll) + E((L, L) T ) (4.12) 



by means of the predictable quadratic variation of the process L, cf. nT3l 1.4.2 and 1.4.50(b)]. 
Using the integral structure of payoff and strategy, we obtain by a stochastic Fubini argu- 
ment that 

pR+ioo pR+ioo 

H-c-cp'S T = S T p(z)dz-c- (<p(z) • St) p(z) dz. (4.13) 

J R—ioo J R—ioo 

The key idea is to identify a family of martingales l(z), z E R + iR, such that 

l(z) T = S z t - cp(z) • S T . 

Then the process 

R+ioo 

l(z)p(z)dz-c (4.14) 

R—ioo 
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is the canonical candidate for the martingale L, and the bilinearity of the predictable covari- 
ation (•, •) suggests that 

/ pR+ioo pR+ioo 

(L,L) = n l(y)p(y)dy, l{z)p{z)dz 

\ J R—ioo J R—ioo I ,a c-\ 

R+ioo pR+ioo ^ ' ' 

I (Kv)>K z )) p{y)p{ z )dydz. 

R—ioo J R—ioo 

In this case, Fubini's Theorem implies 

pR+ioo pR+ioo 

E((L,L) T )= / / E((l(y),l{z)) T )p{y)p{z)dydz. 

J R—ioo J R—ioo 

We now consider how to determine l(z). If the stock price S is a martingale, the definition 
of the cumulant generating function and the martingale property of cp(z) • S yield that 

E (S Z T - <p(z) • S T \&t) = ^e"W< T -*> - <p(z) • St 
is the appropriate martingale. Motivated by this fact, we make the ansatz 

l{z) t = S z t a{z,t)-(p{z)-S t (4.16) 

for deterministic functions a(z) : [0, T] — > C with a(z,T) = 1 in the general case. The 
drift rate of (4.16 ) can be calculated using integration by parts, and setting it to zero yields a 
linear ODE for the mapping t h-> a(z, t). The solution then leads to the desired candidates 
for l(z) and L via (|4.16[) resp. (|4.14|). 



5 Numerical illustration 

In this section, we illustrate our formulas by two examples. First, we compare the perfor- 
mance of the Black-Scholes strategy and the variance-optimal strategy in the normal inverse 
Gaussian (henceforth NIG) Levy model (cf. 0) with parameters inferred from a historical 
time series, i.e., we use the physical probability measure. Afterwards, we assess the hedging 
errors of the Black-Scholes strategy, the delta hedge and the variance-optimal strategy in the 
diffusion- extended CGMY (CGMYe) model (cf. [4j) for parameters obtained by a calibra- 
tion to market prices of options, i.e., in this case we use a risk-neutral probability measure. 
To evaluate the mean- variance optimal hedging error and the corresponding optimal initial 
capital, we use the formulas of [fTTTl . 

We assume a riskless interest rate of 4% and consider a European call option with ma- 
turity T = 0.25 years and discounted strike K = 99. The integral representation of the 
corresponding payoff function is given in Example \2.2\ 

For the numerical evaluation of the integrals we use the routine gsl_integration_qagi 
from the quadrature framework of the GNU Scientific Library (GSL, cf. [9]), which com- 
putes integrals on the whole real line by a transformation on [— 1, 0) U (0, 1] and then applies 
an adaptive algorithm using a 15-point GauB-Kronrod scheme. For the parameter R we 
choose R = 1.1. The outcome of the numerical computations is quite robust with respect to 
the choice of R; instabilities occur only if R is chosen very close to or very far away from 1. 
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1.6 




Initial stock price 
(b) Relative hedging error 



Figure 1 : Hedging errorr for varying initial stock price (NIG model under physical measure) 
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Drift rate 

Figure 2: Mean squared hedging error for varying drift rate /t NIG (l) (NIG model under 
physical measure) 



5.1 Physical measure 

In the following example, we examine the quality of the Black-Scholes strategy as a proxy 
for the variance-optimal hedge in the NIG model. Recall that by Lemma |3.5[ the Black- 
Scholes strategy is given by the replicating Black-Scholes delta using the volatility pa- 
rameter a of the physical price process dynamics, i.e., no risk-neutral (pricing) measure 
is necessary to determine the hedging strategy in this case. 

The cumulant generating function of the NIG Levy process is given by 



(z) = fiz + 5 (V" 2 - /? 2 - V® 2 + z) 2 ) 



for n e M, 5 > 0, < \(3\ < a and z G {y G C : \(3 + Re (y) | < a} (cf., e.g., OH 
Section 5.3.2]). As for parameters, we use 

a = 75.49, /3 = -4.089, 5 = 3.024, fj, = -0.04, 



which corresponds to the annualized daily estimates from 11221 for a historical time series 
of Deutsche Bank, assuming 252 trading days per year. It is easily verified that this market 
satisfies the prerequisites of Section [2] The volatility parameter a for the Black-Scholes 



strategy (cf. Lemma 3.5) is set to 0.2 so that the log-returns in the corresponding Black- 
Scholes market and in the NIG Levy market exhibit the same variance. 

Figure [T] shows the mean squared hedging error and the relative hedging error (i.e., the 
root of the mean squared hedging error divided by the initial capital) of the Black-Scholes 
strategy and the mean-variance optimal strategy for varying initial stock price. As initial 
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capital for the Black-Scholes hedge we use the Black-Scholes price of the option, which is 
virtually indistinguishable from the variance-optimal initial capital, though (compare ifTTTO . 
For an at-the-money call the relative hedging errors differ by 4.73% and amount to 0.113 
(mean- variance optimal) and 0.118 (Black-Scholes). 

Figure [2] illustrates how the two strategies react to different drift rates of the underlying. 
More specifically, the figure shows the mean squared hedging error for an at-the-money call 
option with strike K = 100 and maturity T = 0.25 for varying drift rate «(1), controlled by 
varying the location parameter /i of the NIG process. Since the Black-Scholes strategy does 
not incorporate such systematic drifts directly, the two errors differ least in the martingale 
case k(1) = 0. Altogether, the Black-Scholes strategy seems to be a surprisingly good 
proxy for the mean- variance optimal hedge, particularly for moderate drift rates. 



5.2 Risk-neutral measure 

In this second example, we compare the performance of the Black-Scholes hedge, the delta 
hedge and the variance-optimal hedge in the CGMYe model. The determination of the delta 
hedge and the variance-optimal hedge as well as the computation of the hedging error all 
take place with respect to a risk-neutral measure inferred from a calibration of the model 
to option prices. The volatility parameter a for the Black-Scholes hedge is chosen for each 
initial stock price such that the resulting Black-Scholes price of the option matches the one 
implied by the CGMYe model. 

Let us emphasize that one has to be careful with the interpretation of hedging errors 
computed under a risk-neutral probability measure. This is because expected values under 
this measure have no direct statistical meaning. Nevertheless, this approach is quite com- 
mon in the literature, cf., e.g., [El and the references therein. For a more detailed empirical 
investigation taking care of this issue, one would first determine the delta hedge under a 
risk-neutral measure obtained by calibration. In a second step, one would then estimate a 
parametric ansatz for the market price of risk in order to switch to an appropriate physical 
measure. However, this is beyond our scope here and is therefore left to future research. 

The cumulant generating function of the risk-neutral CGMYe Levy process is given by 

=UJ _ 1^ + 1^2 + C p ( _ y) ((M _ z) Y _ M Y + {Q + g) Y _ Q Y} ? 

where Y denotes the Gamma function and C > 0, G > 0, M > 0, Y < 2 and r] > 
are the parameters of the model. Given these quantities, the value of to is chosen such that 
K CGMYe (l) = 0, which implies that the stock price process is a martingale. For a more 
detailed discussion of the meaning of the different parameters we refer to 01 . Here, we use 
the values from the same paper given by 

C = 9.61, G = 9.97, M = 16.51, Y = 0.1430, 77 = 0.0458. 

Note that since t] > 0, the driving Levy process has a Brownian component and hence 



the delta hedge exists and is a A-strategy by Remark 3/7 Moreover, the model satisfies the 
prerequisites of Section |2j 
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Figure 3: Hedging errors for varying initial stock price (CGMYe model under risk-neutral 
measure) 



14 



Figure [3] shows the mean squared hedging error and the corresponding relative hedg- 
ing error (cf. Section 5.1 ) of the variance-optimal strategy, the Black-Scholes and the delta 
hedge for different initial stock prices in the CGMYe model with parameters as explained 
above. The mean squared errors for an at-the-money option amount to 12.57 for the variance- 
optimal hedge, to 14.68 for the Black-Scholes hedge and to 16.41 for the delta hedge. This 
corresponds to a deviation of the relative errors of 8.10% (Black-Scholes) and 14.3% (delta) 
from the variance-optimal value. 

In contrast to the example in Section |5.1[ mean squared and relative hedging errors in 
general are much higher and the Black-Scholes strategy performs, compared to the variance- 
optimal one, worse than in the previous example. This can be explained by the fact that 
skewness and especially excess kurtosis of the driving Levy process are much more pro- 
nounced in this case. Indeed, in the CGMYe model skewness and excess kurtosis of the 
daily logarithmic returns amount to -3.852 and 62.32 resp. to -0.2384 and 0.2416 for the 
yearly logarithmic returns. In contrast, skewness and excess kurtosis in the NIG model are 
given by -0.1709 and 3.356 for the daily logarithmic returns and by -0.0108 and 0.0133 for 
the yearly returns. Moreover, the mean squared and relative hedging error of the Black- 
Scholes strategy is considerably lower than that of the delta hedge. This effect has also been 
observed by 0J for stochastic volatility models in a discrete time setup. 

Altogether, despite the much higher skewness and kurtosis, the Black-Scholes hedge 
still seems to be a quite stable proxy for variance-optimal one and, in particular, performs 
noticeably better than the model delta. 



6 Conclusion 

In general exponential Levy models, we have derived semi-explicit formulas for the mean- 
squared hedging error of a European- style contingent claim. This has been done for so- 
called A-strategies, which include the Black-Scholes hedging strategy and more general 
delta hedges. Numerical examples obtained by implementing these results show that - both 
under the physical and subject to a risk-neutral probability - the Black-Scholes hedge seems 
to perform surprisingly well also in Levy models with jumps. Nevertheless, it does lead to 
a non-trivial hedging error in this case, which can be quantified using our approach. 



A Proof of the main result 



In this appendix, we present the proof of Theorem 4.2, which is split up into several in- 
termediate statements. An essential tool for the forthcoming considerations are the special 
semimartingale decomposition and the predictable covariation of complex powers of S, 
provided by the following 

Lemma A.l. For z G R + iR, the process S z is a special semimartingale whose canonical 
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decomposition S z = Sq + M(z) + A(z) is given by 

M(z) t = [ e Kiz)s dN(z) s , A(z) t = K(z) I S z s _ds, 



o 



where 

N(z) t := e- K{z)t S z t . 

Moreover, for y, z G R + iM. and continuously differentiable functions (3, 7 : [0, T] — > C, 
the process [S y (3, S 2 ^] is a special semimartingale with compensator (S y (3, S Zr y) given by 

(Sy/3, S z 1 ) t = ( K (y + z)- k{ V ) - k(z)) f Stasis ds. (A.17) 



PROOF. This follows along the lines of the proof of [|TTL Lemma 3.2] . □ 

With the special semimartingale composition of S at hand, we can now establish that 
the mean squared hedging error of a A-strategy is well-defined. 

Lemma A.2. The /^.-strategy <p satisfies ip G L(S) and ip • St G L 2 (P). 

Proof. Fubini's Theorem yields the predictability of cp. The assertion then follows from 
[fTTl Lemma 3.1], Fubini's Theorem, Holder's inequality and Condition 2 of Definition [3TT 



since t ^ E(\S z t | 2 ) = S$ R e tK W is bounded on [0,T] for z G R + zR by fll4] >. □ 

The following proposition ascertains that deterministic integration in the representation 
of cp and stochastic integration with respect to S can be interchanged, compare ( |4.13[ ). 

Proposition A.3. We have 

/•R+ioo 

ip • S = I ((p(z) • S) p(z) dz. 

J R—ioo 



PROOF. By the definition of L(S) in Q3] 111.6.6(c)], Lemma A.l and since E(S 1 



g2 e tK{2) anc j me | oca Uy bounded process £_ are bounded resp. pathwise bounded on [0, T], 
a deterministic process (H t )te[o,T] belongs to L(S) if J Q T Hfdt < 00. Hence Condition 2 of 



Definition |3 . 1 1 implies in combination with Fubini's Theorem that 

-R+ioo \ 1 / 2 



(pR+100 \ 
/ |s(z,.)| 2 |p(*)||^| eL(S). 
JR-ioo J 



Since S^ 1 is locally bounded, 031 111.6.6.19(e)] yields 

-R+ioo \ V 2 / r R+ioo \ 1 / 2 



(/•R+too \ / rR+ioo \ 

/ Hz)\ 2 \p(z)\\dz\) =S R -U \g(z,.)\ 2 \p(z)\\dz\) eL(S). 

JR-ioo / \JR-ioo J 



The assertion now follows from Fubini's Theorem for stochastic integrals, cf., e.g., Il20l 
Theorems 63 and 65]. □ 
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The following theorem shows that our Ansatz ( |4.16[ ) indeed works. 
Theorem A.4. Let z G R + iR and define a(z) : [0, T] — > C as 

a{z,t):= (l-K(l) C ' e K ^ s ~ T) g(z, s) ds] e "W( T -*). 



(A.18) 



(A.19) 



(A.20) 



Then a solves the terminal value problem 

a'(z, i) + k(z)oc(z, t) - K(l)g(z, t) = 0, t G [0, T], 
«0,T) = 1, 

and ?/ze process l(z) defined by 

l{z) t := S"a(z,t) - (p(z) • S t 
is a local martingale with 1(z)t = — <f(z) • St- 

Proof. Let z E R + iR. Differentiation shows that a solves (A.19). To prove the second 
part of the assertion, we decompose l(z) into a local martingale and a drift, and then con- 
clude that the latter vanishes due to the choice of a. Integration by parts and [fT3l 1.4.49(d)] 
yield 

l( z ) = S z a(z,0) + I S z s _a'(z,s)ds + a(z,-) • S z - <p(z) • S. 
Jo 

It now follows from Lemma [Aj] that 

l(z) = S£a(z,0) + / S*_a'(z, s)ds + a(z, ■) • A(z) + a(z, •) • M(z) 
Jo 

-<f(z) 'A(l)-<f(z) • M(l) 
= S$a(z, 0) + a(z, •) • M(z) - <p(z) • M(l) 



(A.21) 



+ / S z _ (ot'{z, s) + n(z)a(z, s) — n(l)g(z, s) j ds . 



Since a satisfies (A.19), the last integral on the right-hand side of (A.21 ) vanishes. The 
remaining terms are local martingales by 021 1.4.34(b)], because oc{z, ■) and ip(z) are locally 
bounded. Since a(z, T) = 1, this proves the second part of the assertion. □ 



Lemma A.l now allows us to compute the predictable quadratic covariations (l(y), l(z)). 

Proposition A.5. For all y,z G R + iR, the process [l(y), l(z)] is a special semimartingale 
with compensator (l(y), l(z)) given by 

(l(y),l(z)) t = (n(y + z) - n{y) - k(z)) / S y s + Z a(y, s)a(z, s) ds 



- (k{v + 1) - n(y) - k(1)) / S y s + Z a(y,s)g(z,s)ds 

Jo 

- (k(z + 1) - k(z) - k(1)) / S y s + Z a(z, s)g(y, s) ds 

Jo 

+ («(2) - 2«(1)) / S y + Z g(y,s)g(z,s)ds. 



(A.22) 
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PROOF. Let y,z e R+iR. By the definition of /(•) in (A. 20 ), the bilinearity of the quadratic 
covariation [•, •] and flT3l 1.4.54], we obtain that 

[l(y),l(z)} = [Sya(y) } S z a(z)}- V (z) • [5»a(y),5] 

-cp(y).[S,S z a(z)) + (ip(y)ip(z))-[S,S], 



because ip(y) and <p(z) are locally bounded. Recall that by Lemma A. 1[ the square bracket 



processes on the right-hand side are special semimartingales with compensators given by 
( |A.17[ ). Again using that ip{z) is locally bounded, it then follows from [[131 1.4.34(b)] that 

[l{y),l(*)] ~ (S y a(y),S z a(z)}- V (z) • (S y a(y) 7 S) 

- V (y)-(S,S z a(z)) + ^(yMz))-(S,S) 

is a local martingale. By inserting the explicit representations ( |A.17] >, we obtain that our 
candidate for (l(y), l(z)) indeed compensates [l(y), l{z)]. Since it is also predictable and of 
finite variation, this completes the proof. □ 

The following technical lemma provides the p-integrability for several expressions, which 
is necessary to apply Fubini arguments in the proofs of Propositions A. 8 and |A.9 



Lemma A.6. With a and l(z) defined as in (A.18) and (A.20), we have the following for all 
te[0,T}: 

1. There exists a constant b\ such that Re (k(z)) < b\for all z G R + iR. 

2. The mappings z t— > a(z, t) and z y a(z, t) 2 are p-integrable. 

3. The mapping (y, z) \-t E (\(l(y), l(z)) t \) is twice p-integrable. 

4. The mapping (y, z) \-t (l(y),l(z)) t (u) is twice p-integrable for almost all u e Q. 
Proof. Let y, z e R + iR and t E [0, T). 

1. The definition of the cumulant generating function and Jensen's inequality yield 



3 Re(K( 2 )) 



\E (e 2Xl )| < E (|e 2Xl |) = E ( e Re ^ Xl ) = e K 



(Re(z)) 



< e 



In 



forfex := k(R) V0. 
2. By Holder's inequality, we have 

Hz,t)\ <e Re ^ T ^ + \K(l) 



e Re( K(2 ))( S -t) |^ (ZjS) 



ds 



< e 6lT + |«(1)| e^-'^T - t)* ^ \g{z, s)\ 2 ds^j 



1/2 



The p-integrability of z (->■ a(z,t) now follows, because z !->• 1 ist p-integrable by 



Assumption |2. 1| and the last integral on the right-hand side is p-integrable by Holder's 
inequality and the fact that Condition [2] in Definition 3.1 holds for the A-strategy if. 



Considering the square of the right-hand side, it follows directly from Assumption 2.1 



and ConditionSin Definition 3.1 that z \-t a(z, t) 2 is p-integrable 
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For the proof of Assertions [3] and [4j first note that the bilinearity of (•, •) yields 

±Re «%), l(z))) = \ ( (%) ± W), Ky)±W)) ~ (%),%)} - (Kz),W)) ) (A-23) 
and 



l(y) ± l(z), l(y) ± l(z) ) < (l(y) ± l(z), l{y) ± l(z)) + (l(y) T l(z),l(y) =F /(*) 



Applying an analogous polarization argument to lm((l(y),l(z))) by replacing l(z) with 



il(z), we see that it is sufficient to consider only the covariation of the form \ l(z), l(z 
in order to show [3] and |4j Since this process is real-valued and increasing, we can restrict 
ourselves to t = T. 



3. By applying the arguments from the proof of Proposition A. 5 
[l(y), l{z)\ and taking absolute values, we obtain 

f-T 



to 



l(z),l(z) 



instead of 



(l(z),l(z)) T < J S™ c{z) \a(z,s)\ 2 \K(2Re(z))-2Re(K(z))\ ds 
+ [ T S* Rc(z) \g(z,s)\ 2 \K(2)-2K(l)\ds 



Using 



+ / 2S 2 ^ (z) \a(z,s)\ \g{z,s)\ \k(z + 1) - k{z) -k(1)| ds. 
Jo 

(A.24) 



E (s 2Rc{z) ) = E (Sf*M) = 5 JR«W e «(2B«W) < S 2R ( e T K (2 R) v ^ 



and Fubini's Theorem, we obtain from ( |A.24| ) that 



E^l(z),l(z)^J < (5 2fl (e TK(2ii) Vl)) (J \a(z,s)\ 2 \K(2Re(z)) -2Re (k(z))\ ds 



+ / \g(z,s)\ 2 \ K (2) -2k(1)| 



+ y 2|a(*,s)| |#0,s)| |k(z + 1) - «(*) - «(1)| dsj . 

(A.25) 

To prove Assertion 3, it therefore suffices to show that all integrals in ( |A.25[ ) are p- 
integrable. For the first one, we have 



|aO,s)| 2 K2Re(V)) -2Bb(k(z))\ ds 

r T r T 

<\k{2R)\ \a{z,s)\ 2 ds + 2 / \a(z, s)\ 2 |Re {k{z))\ ds (A.26) 
Jo Jo 
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Let us first consider the second integral in ( |A.26[ ). Inserting the representation ( |A.18[ ) for 
a, we obtain 



< 



| Re (k(z))\ \ce(z, s)| 2 ds 

|Re(/t(z))|e 2Rc(K(z))(T - s) ds 

|Re(/t(z))|e Rc(K(2))(T - s) |k(1)| 

eT 



3 Re(«( Z ))(r-s) 



g{z, t)| dr J ds 



+ 



< 



+ 2 



|Re(K(z))| |K(1)| 2 N e ReW,))(r-s) 
Re(;2) | e 2ReW,))(T- S ) ^ 

|Re(/€(z))|e Re(K(2))(T - s) |fi:(l)| 



^(z, r)| dr) ds 



+ 



r 



\g(z, r)| dr 1 ds 



|Re(ic(z))||/c(l)| s 



,2Re(«(z))(r-8) rfr 



^ |fi-(2,r)| 2 dr^j ds, 



where we applied Holder's inequality twice in the last step. Using that Re (k(z)), z £ 
R + iR, is bounded from above by the constant b% > 0, it is easily seen by elementary 
integration that the integrals of the form 

/ \Re(K(z))\e mRe{K{z))s ds, < a < T, me {1,2}, 

J a 

are uniformly bounded on R + iR. Moreover, the terms \g(z, r)| 2 dr are bounded 
by J„ T \g(z, r)| 2 dr, which is p-integrable since (p is a A-strategy. The p-integrability of 



J | r) | dr J follows after an application of Holder's inequality. Altogether, this 
yields that the second integral in ( |A.26[ ) is p-integrable. The p-integrability of the first 
integral in ( |A.26[ ) follows analogously by inserting the representation for a from ( |A.18| ) 
and exploiting that <p is a A-strategy. By the latter fact, we obtain directly that the second 
integral in ( |A.25| ) is p-integrable as well. To deal with the third one, we use the inequality 



\k(z + 1) - k(z) - k(1)| 2 < («;(2Re (z)) - 2Re («(z))) («(2) - 2k(1)) 
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established in ifTTl Lemma 3.4] and apply Holder's inequality to conclude that 



\a(z, s)\ \g(z, s)\ \k(z + 1) — k(z) — ds 



< 



T 

\a(z, s) 



ds 



|/c(2Re (z)) - 2Re (k(z))\ \k(2) - 2«(1)| ds 



o 



Using Holder's inequality, this proves Assertion 3, because the squares of both integrals 
on the right-hand side have already been shown to be p-integrable. 

4. Note that in ( |A.24[ ) the only stochastic terms are given by S^^ z ' = S 211 . Since this 
process is locally bounded, almost all of its paths are bounded on [0, T]. The estimates 
from the proof of Assertion [3] therefore also show that As sertion|4] holds. □ 

The estimates in Lemma [A3 immediately yield the following 



Corollary A.7. For all z e R + iR, the process l(z) defined in (A.20) is a square-integrable 
martingale. 



PROOF. Recall that l(z) has already shown to be a local martingale in Theorem A.4 Since 
[l(z), l(z)] — (l(z), l(z)) is a local martingale as well, the assertion follows from [[131 1.4.50(c)] 
by localization, monotone convergence and the fact that the right-hand side of ( |A.25[ ) is fi- 
nite for all z e R + iR. □ 



The next two propositions show that the candidates proposed in Equations ( |4.14[ ) and ( |4.15[ ) 
indeed coincide with the desired martingale and its quadratic variation of Ansatz ( |4.12[ ). 

Proposition A.8. The process L defined by 



L, :-- 



l(z)tp(z) dz — c 



(A.27) 



R—ioo 



is a real-valued square-integrable martingale with Lt = H 



ip • St- 



Proof. First note that by Lemma A. 6 2 ) and Proposition A. 3 the integral in (A.27 ) is well- 
defined. Fubini's Theorem and dominated convergence show that f^-™ S^a(z } t) p(z) dz 
is an adapted cadlag process. For t G [0, T\, Holder's inequality and another application of 
Fubini's Theorem yield 



E 



R+ioo 



l(z) t p(z) dz 



R—ioo 



< E 



R.+ioo 



\l(z) t \ 2 \p{z)\ \dz\ 



R—ioo 
R+ioo 



R—ioo 



E{\l(z) t \ 2 ) \p(z)\ \dz\ 



R+ioo 

\p(z)\ \dz\ 

R—ioo 

R+ioo 



\p(z)\ \dz\ 



R—ioo 
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a(z,0)\ 2 + E^l(z),l(z)^ T y (A.28) 



Since l(z) is a square-integrable martingale by Corollary A.7| |fT3l 1.4.2 and 1.4.50(b)] imply 

E(\l(z) t \ 2 ) = \l(z) \ 2 + E^l(z),T^))^ <£*«'-■'- - |2 

Because the right-hand side of ( |A.28[ ) is p-integrable by Lemma A.6[2p ) and independent 
of t G [0, T\, we conclude that L is indeed square-integrable, i.e., 

sup E (|£t| 2 ) < oo. 

te[o,T] 

In order to show the martingale property of L, consider arbitrary < s < t < T and 
F G =f s . By Fubini's Theorem and the martingale property of l(z), we have 

-R+ioc 



E((L t -L s )l F ) = E 



(l(z) t - l(z) s ) l F p{z) dz 

R—ioo 
R+ioo 

E((l(z) t -l(z) s )l F ) p(z) dz = 0, 

R—ioo 



and hence E(L t \J^ s ) = L s . Since a(z, T) = 1, it follows from Assumption 2.1 and Propo- 



sition 



A. 3 that Lt is given by the asserted, real-valued random variable. The martingale 



property of L then yields that L t is real- valued for all t G [0,T], which completes the 
proof. □ 



Proposition A.9. The predictable quadratic variation (L, L) of the process L defined in (A.27) 
is given by 



(L,L) t 



R+ioo pR+ioo 



R—ioo J R—ioo 



(Ky)A z ))t p(v)p( z ) d v dz - 



(A.29) 



PROOF. By Proposition A. 8 and [13, 1.4.2 and 1.4.50(b)], it suffices to prove that the can- 



didate C := fn^i™ Ir- 1 ^ (Kv)i K z )) p(u)p( z ) dydz is a predictable process of finite vari- 
ation such that I? — C is a local martingale. First, note that the integral in ( |A.29[ ) is well- 



defined by Lemma |A^6j|4]). Moreover, since t i-> (l(y),l(z)) t is continuous by Proposi- 
tion|A.5|and hence predictable for all y, z G R + iR, the process C is predictable as well by 



Fubini's Theorem. To see that C is of finite variation, note that it is a linear combination of 
expressions of the form 

-R+ioo /•R+ioo 

Re{(l(y),l{z))) h(y)j(z)dydz 



pH+ioo pi 
J R-ioo J R- 



R—ioo J R—ioo 
R+ioo pR+ioo 



or 

/ / lm((l(y),l(z)))h(y)j(z)dydz 

J R—ioo J R—ioo 

for h, j G {Re (p) + , Re (p)~ , Im (p) + , Im } . In view of ( |A.23[ ), we obtain that 

-R+ioo pR+ioo 

Re({l(y),l(z)}) h(y)j(z)dydz 



R—ioo J R—ioo 

pR+ioo /-R+ioo 

2 



l(y) + l(z), l(y) + l(z) ) h(y)j(z) dydz 



R—ioo J R—ioo 

R+ioo f R+ioo 



R—ioo J R—ioo 



l(y),l(y)) + (l(z),l(z))) h(y)j(z)dydz 
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which is the difference of two increasing adapted processes and hence of finite variation. 
The argument applies analogously to lm((l(y),l(z))). To show the martingale property 
of L 2 — C, we can assume without loss of generality that c = 0. Observe that 

-R+ioo pR+ioo 

(l{y) t l{z) t -(l(y),l{z)) t )p(y)p(z)dydz 



(pR+ioo p 
/ / 
J R—ioo J R 



R—ioo 
R+ioo pR+ioo 



J R—ioo J R—ioo 

Moreover, E(\l(y) t l(z) t \) < E{\l{y) t \ 2 ) l l 2 E{\l{z) t \ 2 fl 2 and E(\(l(y),l(z)) t \) are twice 
p-integrable by Lemma |A.6p j), Holder's inequality and the fact that the right-hand side 
of ( |AT28] ) is p-integrable. This shows that L 2 - C t E L\P). For < s < t < T and 
F E J^s, we can therefore apply Fubini's Theorem to obtain 

pR+ioo pR+ioo 

E ((L 2 - Ct) 1 F )= / E {{l{y)tl{z)t - (Kv),K*))t) lr) P(y)p(z)dydz 

J R—ioo J R—ioo 



R+ioo fR+ioo 



E ((l(y) s l{z) s - (l(y), l{z)) s ) 1 F ) p(y)p(z) dydz 



R—ioo J R—ioo 

E((L 2 S -C S ) 1 F ) 



where we used in the second step that l(y)l(z) — (l(y), l(z)} is a martingale by Corollary A. 7 
and [SI 1-4.2 and 1.4.50(b)]. Hence, E (L 2 - C t \ = L 2 S - C s , which completes the 
proof. □ 

Now we can finally prove our main Theorem |4.2| by combining the preceding results. 

Proof of Theorem 14 .21 In view of Lemma \A72\ the mean squared hedging error corre- 
sponding to the endowment/strategy pair (c, cp) is well-defined. By Assumption 2.1 Propo- 
A.3 the definition of l(z) in ( |A.20[ ) and Proposition |A.8[ we have 



sition 



E((H -c-<p- S T ) 2 ) = E(L 2 T ) = E(L 2 ) + E((L, L) 



where the second equality follows from Hi 31 1.4.2 and 1.4.50(b)]. Now notice that by defini- 
tion, 



r2 



R+ioo "\ 2 

SqO.(z, 0) p(z) dz — c 



R—ioo 



In view of Lemma A. 6 3), Fubini's Theorem and Proposition A. 9 yield 

-R+ioo /•R+ioo 



E((L,L), 



R—ioo J R—ioo 



K z )) T )p(y)p( z ) d v dz - 



Since E (S^ z ) = S^ +z e tK<yV+z \ and because the continuous functions t \g{z,t)\, t i-> 
|a(z,i)| and 

t^E(\S£*\) <S 2R e tK ^ 



are bounded on [0, T], Proposition A. 5 and another application of Fubini's Theorem com- 
plete the proof. □ 
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